1 



HIGH ENERGY EMISSION FROM IGR J16320-4751 



L. Foschini 1 , J. A. Tomsick 2 , J. Rodriguez 3 \ R. Walter 4 , A. Goldwurm 3 , S. Corbel 3 5 , and P. Kaaret 6 

^ASF-CNR, Sezione di Bologna, Via Gobetti 101, 40129 Bologna (Italy) 
2 Center for Astrophysics and Space Sciences, University of California at San Diego, La Jolla, CA 92093, USA 
3 CEA Saclay, DSM/DAPNIA/SAp (CNRS FRE 2591), 91191, Gif-sur-Yvette Cedex, France 
4 INTEGRAL Science Data Centre, Chemin d'Ecogia 16, 1290 Versoix, Switzerland 
5 Universite de Paris VII, Federation APC, 2 place Jussieu, 75251 Paris Cedex 05, France 
6 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 



ABSTRACT 



IGR J16320 - 4751 was re-discovered by IBIS/ISGRI 
on board INTEGRAL in early February 2003 during the 
observation of the black hole candidate 4U1630 — 47 (PI 
Tomsick). This source, already observed by ASCA and 
BeppoSAX, belongs to the class of heavily absorbed ob- 
jects (Nh > 10 23 cm -2 ) that populate some arms of the 
Galaxy. Soon after the rediscovery by IBIS/ISGRI, the 
source was observed by XMM-Newton: the arcsec po- 
sition found with XMM allowed ones to find the most 
likely infrared counterpart. We present here the reanaly- 
sis of the high energy emission from IGR J16320 — 4751 
detected by IBIS/ISGRI, including the spectral and tem- 
poral characteristics. We also present a reanalysis of the 
XMM-Newton and optical/IR data. 

Key words: X-rays: binaries; X-rays: individuals: IGR 
J16320-4751. 



1. INTRODUCTION 



IGR J16320 — 4751 was serendipitously discovered on 
February 1.4, 2003 UT, with the IBIS/ISGRI detector 
(Ubertini et al. 2003, Lebrun et al. 2003) on board 
the INTEGRAL satellite during the AOl observation of 
the black hole candidate 4U 1630 - 47 (PI Tomsick). 
The coordinates (J2000) were a = 16 h : 32 m .0 and 
S = —47° : 51 m .0, with an uncertainty of 2' (Tomsick et 
al. 2003a). A quick archival research allowed ones to find 
immediately the X-ray counterpart in the ASCA catalog, 
as the source AX J1631.9 - 4752 (Sugizaki et al. 2001). 
IGR J16320 — 4751 is therefore apparently located in the 
Norma Arm of the Galaxy, an active star forming region 
distant from the Sun about 5 kpc (Georgelin & Georgelin 
1976, Russeil 2003). 

A Target of Opportunity request (ToO) with XMM- 
Newton was immediately activated and it found one sin- 
gle source inside the error circle of ISGRI. The coor- 
dinates (J2000) were a = 16 h : 32 m : 01 s .9 and 



5 = -47° : 52 m : 29 s with an uncertainty of 4" (Ro- 
driguez et al. 2003a). In this reduced error circle it was 
possible to identify two reliable infrared counterparts in 
the 2MASS catalog (Tomsick et al. 2003b). 

A preliminary analysis of the XMM-Newton observation, 
together with the identification of the infrared counter- 
part was presented in Rodriguez et al. (2003b). Here we 
present a reanalysis of the INTEGRAL AOl observation, 
together with a reanalysis of the XMM-Newton ToO data. 



2. INTEGRAL DATA ANALYSIS 



INTEGRAL observed the Norma region from 1 February 
2003 05 h : 40 m : 58 s to 5 February 2003 07 h : 53 m : 01 s 
UT for an elapsed time of 300 ks and with the dither pat- 
tern 5x5. The screening, reduction, and analysis of the 
IBIS/ISGRI data have been performed by using the IN- 
TEGRAL Offline Scientific Analysis (OSA) v. 3 (Gold- 
wurm et al. 2003a), available to the public through the 
INTEGRAL Science Data Centre 1 (ISDC). 



2.1. Imaging 

An initial analysis was performed by deconvolving the 
ISGRI shadowgrams to obtain images. IGR J16320 — 
4751 was detected in the imaging pipeline in the energy 
bands 20 - 40 keV (signal-to-noise ratio SNR = 30a) 
and 40 — 60 keV (SNR = 11a). The mosaic images are 
shown in Fig. [l] Given the present uncertainties in the 
validation of the software for the off-axis sources, it was 
adopted the same procedure described by Goldoni et al. 
(2003). We divided the count rate of IGR J16320 - 4751 
by the count rate of the Crab in a similar off-axis angle, 
extracted from the calibration observations. We added a 
5% of systematic error, to take into account the residual 
fluctuations in the count rate (cf Goldwurm et al. 2003b). 
The calculated flux is (8.0 ± 0.5) x 10" 11 erg cm -2 s _1 

^ttp : //isdc . unige . ch/ index . cgi?Sof t+download 
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Figure 1. IBIS/ISGRI mosaic of the Norma region centered on the BH 4U1630 — 47. (left) Energy band 20 — 40 keV, 
(right) Energy band 40 — 60 keV, with a total exposure of 297 ks. Equatorial coordinates are superimposed. 



and (2.0±0.2) x 10 -11 erg cm -2 s _1 , in the energy bands 
20 — 40 keV and 40 — 60 keV respectively. No significant 
detection was recorded at higher energy. 



2.2. Timing 



IGR Jl 6320— 4751 was detected in most of the individual 
ScWs with 3 < SNR < 6a. Therefore, to study the vari- 
ability of the source it was decided to rebin the lightcurve 
so that each time bin contains 8 ScW, corresponding ap- 
proximately to 15 — 20 ks. The obtained lightcurve shows 
a clear outburst at the time of the discovery, plus some 
other periods with a certain variability (Fig. |2|>. During 
the outburst, the flux in the band 20 — 40 keV raised from 
about 7-8 mCrab to 30 mCrab (Fig.|3 left). The time 
scale of the variations occurs on w 10 4 s or even more. 
A similar behaviour is clearly seen in the 40 — 60 keV 
energy band (Fig. right). 

Since in a coded-mask instrument the brightest sources 
in the field of view (FOV) can - under certain conditions 
- significantly affect the detection of the other nearby 
sources, we investigated also the time variability of the 
three brightest sources in the FOV, namely 4U1630 — 47, 
4U1700 - 377, and GX340 + 0. All the analysed sources 
display different variability patterns and therefore we 
conclude that the variability of IGR J16320 - 4751 is 
genuine. 



2.3. Spectral extraction 

To extract the spectrum of the whole observation and also 
during the outburst we used both the spectral extraction 
pipeline and the results from the image analysis, to per- 
form a check of the results. This is necessary to take into 
account that IGR J16320 — 4751 is located in a crowded 
region, with at least two nearby sources (4U 1630 — 47 
and IGR J16318 - 4848), and therefore the different ap- 
proach in the deconvolution of images and extraction of 
spectra (see Goldwurm et al. 2003a for a full explana- 
tion of the algorithms) could generate results not always 
consistent each others. 

For the spectral extraction pipeline we rebinned the latest 
RMF matrix from the original 2048 channels to 21 chan- 
nels, by putting all the channels above 200 keV into the 
last channel of the rebinned matrix. The remaining chan- 
nels in the range 15 — 200 keV were grouped in bin of 
about 10 keV in size. This approach has been selected to 
emphasize the energy range 15 — 200 keV and to have 
enough statistics to perform the \ 2 test m m e spectral fit 
with xspec (v 1 1.3). Nonetheless, the source is faint and 
does not allow a fit with multiple component models. We 
decided to fit with a simple power law model to measure 
if there is hardening/softening during the outburst. 

The spectrum of the whole observation has V = 3.8 ±0.9, 
with a flux of about 9 x 10 -11 erg cm -2 s _1 in the energy 
band 20 — 60 keV, consistent with the results from the 
imaging analysis. Also in the spectral extraction pipeline 
no flux is detected for energies greater than w 60 keV. 

The spectrum extracted during of the outburst (bins 2-4 
in Fig. |2> has the photon index varying from Y w 3.1 to 
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Figure 2. IBIS/ISGRI lightcurves oflGR /16320 - 4751 in the energy band 20 - 40 keV (left) and 40 - 60 keV fright/ 
Each bin is composed of 8 ScW, corresponding to about 15 — 20 ks. Upper limits are at 3a level. 




Figure 3. EPIC-PN image in the 0.2 — 12 keV energy 
band, with 5 ks of exposure (i.e. cleaned from soft-proton 
flares). The IBIS/ISGRI error circle is superimposed. It 
is clearly visible only one source inside the error circle. 

T w 2.6, therefore slightly harder than that of the whole 
observation, with a peak flux of 2 x 10 -10 erg cm -2 s _1 
in the 20 — 60 keV energy band. 

During the "normal" activity, investigated only with the 
imaging pipeline, IGR J16320 — 4751 is barely detected 
in the energy band 40 — 60 keV (Fig. O right), with the 
exception of the outburst. 

The hardness ratio (not shown), calculated by using the 
fluxes in the two energy bands (20 — 40 keV and 40 — 
60 keV), is in agreement with the results of the spec- 
tral extraction pipeline, and confirms the softness of the 
source during almost the whole observation. 



3. XMM - NEWTON DATA ANALYSIS 



XMM-Newton observed IGR J16320 - 4751 as a Target 
of Opportunity (ToO) from 4 March 2003 20 h : 20 m : 



38 s UT to 5 March 2003 03 h : 47 m : 18 s UT, with an 
elapsed time of about 26 ks. Due to soft-proton flares, 
the effective exposure was about 5 ks, but sufficient to 
clearly identify (uncertainty < 4") the X-ray counterpart 
at coordinates (J2000) a = 16 h : 32 m : 01 s .9 and 6 = 
-47° : 52 m : 29 s (Rodriguez et al, 2003a,b; see Fig. [3. 
An early analysis of the XMM-Newton data (ToO) were 
published in Rodriguez et al. (2003b). The lightcurve 
in the 5 ks clean from soft-proton flares showed a time 
variability on time scales of 100 s, while the spectrum 
was found to be harder than that of ISGRI, fitted with an 
absorbed power law with photon index T = 1.6^0.1 an d 
N H = (2.lto.t) x 10 23 cm" 2 (x 2 = 55.9, d.o.f.= 69). 
The unabsorbed flux in the 2 — 10 keV energy band was 
1.7 x 10 -11 erg cm -2 s _1 . We refer the reader to the 
paper by Rodriguez et al. (2003b) for more details. 

For the present work, we decided to reanalyze the XMM- 
Newton data to extract the better spectrum. We used still 
the XMM SAS v. 5.4.1 software to process and screening 
the data and the same procedures described in Rodriguez 
et al. (2003b). We adopted a different procedure only to 
extract the spectrum. Since the source is bright enough 
to be clearly visible still in the flared image, we avoided 
the selection of a time region not affected by soft-proton 
flares, and we extracted directly from the contaminated 
image the source plus background counts from a circu- 
lar region centered on IGR J16320 — 4751 with radius 
of 30". Therefore, we extracted the background from 
another circular region with radius 2', and close to the 
source. We performed the background (now including 
the soft-proton flares) subtraction directly in x spec. 

We succeeded to save 17 ks from the elapsed time of 26 
ks, with a loss of only 9 ks (to be compared with the loss 
of 21 ks of the early analysis). We found that the best fit 
model is still an absorbed power law with T = 1.4+q' 2 
and Nh = (2.2 ± 0.2) x 10 23 cm" 2 ( X 2 = 204.5, d.o.f.= 
184), consistent with the results of the early analysis (see 
Fig. |4j. The unabsorbed flux in the 0.2 — 12 keV energy 
band is 2.8 x 10 -11 erg cm -2 s _1 . 

We note that there are some interesting features between 
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Figure 4. EPIC-PN spectrum in the 3 — 12 keV energy band, with 17 ks of exposure and fitted with a single power law 
model. The deviations of data versus model are shown in the bottom window. 



4 and 10 keV, and particularly some hints for emission 
lines of the iron complex, between 6 and 7 keV. The 
addition of a thermal plasma model (mekal model in 
xspec) with a temperature kT — 5.5 keV is able to 
reproduce some of these features, and particularly the 
emission lines of the iron complex. However, this model 
is statistically required only at 91% confidence level. A 
simple gaussian emission line dX E = 6.5 ± 0.1 keV and 
width a = 0.3^0.1 keV is required at 99.67%. Further 
observations are required to better constrain the nature of 
this excess. 



4. OPTICAL/INFRARED COUNTERPARTS 



With the improved position given by XMM-Newton, it be- 
came possible to look for the counterparts at other wave- 
lengths, specifically in the optical/infrared bands (Tom- 
sick et al. 2003b). Two sources have been found in 
the Two Microns All-Sky Survey (2MASS): the first, 
labelled 1 in Fig. |5] could also be the most probable 
counterpart. The infrared colours along the line of sight 
(E(J - H) = 1.0, E(H - K) = 0.8, for the extinction 
due to the Galactic absorption Ay = 11.1) suggesting 
an excess perhaps due to the presence of hot plasma or 
circumstellar dust, consistent with the findings of the X- 
ray analysis. Specifically, the source 1 is not detected in 
the J band, with an upper limit of J > 14. Concern- 
ing the source 2, there are detections in other catalogs: 
a summary of the optical/infrared detections is shown in 
Table[l] It is worth noting that the source 2 is classified as 
a non-star object, that may be either galaxies or blended 
objects, in the Guide Star Catalog 2.2 2 . Indeed, it ap- 

2 http : //www-gsss . stsci . edu/gsc/GSChome . htm 
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Figure 6. Broad band unfolded spectrum of 
IGR /16320 — 4751 obtained with a simultaneous 
fit of XMM-Newton (EPIC-PN) and INTEGRAL 
(IBIS/ISGRI) data. The energy ranges are 3 — 12 keV for 
XMM-Newton and 30 - 60 keV for INTEGRAL. 



pears slightly elongated with an eccentricity of 0.07 and 
a semimajor axis of 3.12 pixels. 



5. DISCUSSION AND CONCLUSIONS 



In the present work we confirmed and extended the re- 
sults obtained by Rodriguez et al. (2003b). The nature 
of IGR J16320 — 4751 is not clear yet and, even though 
the most probable hypothesis is a Galactic X-ray binary, 
the possibility of an extragalactic object cannot be com- 
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pletely ruled out. The source appears to be intrinsically 
absorbed, with an Nh one order of magnitude greater 
than the Galactic absorption along the line of sight. Also 
the infrared observation suggests an excess due to cir- 
cumstellar plasma or dust along the line of sight. The X- 
ray spectrum could support the hot plasma solution, but 
the confidence level is not enough (only 91%) and should 
be confirmed with further observations with much more 
statistics. Therefore, the picture of a high-mass X-ray bi- 
nary (HMXRB) with a compact object accreting through 
winds appears to be very likely for IGR J16320 — 4751. 
The winds could also be the source of obscuration, wrap- 
ping up the compact object. 



The nature of the accreting object is not clear as well: 
the strong differences in the photon index in the energy 
bands 0.5 - 10 keV and 20 - 60 keV could be due to 
changes in spectral states, thus suggesting a black hole. 
On the other hand, the drop of emission for energies 
greater than 60 keV is common in neutron stars. Al- 
though, the INTEGRAL and XMM-Newton are not simul- 
taneous, we tried a simultaneous fit (Fig. 0. The best 
fit model (x 2 = 344.8, dof = 344) is an absorbed 
power law with exponential cutoff (cutoffpl model 
in xspec), typical of accreting neutron stars (see, e.g. 
White et al. 1983). The photon index is T = 0.7±g;|, the 
column density N H = (2.0 ± 0.2) x 10 23 cm" 2 , and the 
cutoff energy 10 ± 3 keV. The scaling constant applied 
to the ISGRI spectrum is 5^. We stress that the obser- 
vations were not simultaneous and the source is strongly 
variable. Therefore, this fit should be taken with extreme 
care. 



Another possibility, given the low luminosity in both 
bands (8 x 10 34 erg/s in the 0.2 - 12 keV energy band, 
and 3 x 10 35 erg/s in the 20 - 60 keV band) if the source is 
located in the Norma Arm (5 kpc), is that we are observ- 
ing the emission from a jet. In this case, we could make 
the hypothesis that the change in the photon index could 
be due to a break, because of a change in the mechanism 
of cooling of electrons. 



The key question in the study of this source is therefore 
the spectral variability. An approved coordinated INTE- 
GRAL and XMM-Newton observation to be performed by 
next August should give us sufficient data to disentangle 
these hypotheses. 
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USNO CCD AC d 


UCAC = 16.0 ± 0.3 
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i?i = 14.6 ±0.3 


USNOBl b 


R 2 = 15.4 ±0.3 


USNO A2 b 


R= 15.0 ±0.3 


GSC 2.2 C 


R= 15.4 ±0.4 
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DENIS e 
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2MASS a 


J = 12.13 ±0.02 


DENIS e 
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DENIS e 
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